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Abstract. - We introduce a model of percolation induced by disorder, where an initially homoge- 
neous network with links of equal weight is disordered by the introduction of heterogeneous weights 
for the links. We consider a pair of nodes i and j to be mutually reachable when the ratio aij 
of length of the optimal path between them before and after the introduction of disorder does 
not increase beyond a tolerance ratio r. These conditions reflect practical limitations of reachabil- 
ity better than the usual percolation model, which entirely disregards path length when defining 
connectivity and, therefore, communication. We find that this model leads to a first order phase 
transition in both 2-dimensional lattices and in Erdos-Renyi networks, and in the case of the latter, 
the size of the discontinuity implies that the transition is effectively catastrophic, with almost all 
system pairs undergoing the change from reachable to unreachable. Using the theory of optimal 
path lengths under disorder, we are able to predict the percolation threshold. For real networks 
subject to changes while in operation, this model should perform better in predicting functional 
limits than current percolation models. 



Percolation theory is the most well-established approaches to study system connectivity, 
and how this connectivity becomes compromised with local system failure [T|. Systems 
usually refer to connected structures such as lattices or random networks [2], and failures 
to the removal of nodes or links. Percolation predicts failure thresholds and size of the 
connected parts of the network after those failures. Practical domains of applicationcs 
include epidemiology [2H5], communication networks like the Internet [B], and propagation 
of information in social networks [7]. 

Recently [S], it was pointed out that there are many contexts in which structural con- 
nectivity may not be sufficient for a network to maintain its functionality when faced with 
failures. Network nodes may need to be more than connected; they may need to be reachable 
from one another, which is practical definition. In regular percolation, node j is reachable 
from node i if there is a path of consecutive links from i to j, i.e., if the nodes are structurally 
connected after any failures. However, in real networks, two nodes initially reachable from 
one another through a structural path can become unreachable after link failures even if 
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there is a new structural path connecting them, but which is in some way unsatisfactory. 
This reasoning lead to the introduction of Limited Path Percolation (LPP), in which reach- 
ability between nodes i and j after the onset of failures is a relative statement: if the path 
length between nodes before (1^) and after {£.■) failures is such that l^/lij < t, then node 
j is reachable from i (and vice versa), where r is an externally imposed tolerance that re- 
flects properties of the functional aspect of the network (regular percolation corresponds to 
r — > oo). For a large class of networks, Ref. [8] shows that a new phase transition appears, 
where networks are now more fragil due to the additional length constraint. 

The approach introduced in Ref. [5] was restricted to the removal (absolute failure) of 
some links, while other links remained unperturbed. However, in many real- world examples, 
links usually do not fail entirely, but can instead become more difficult to use such as in road 
networks affected by weather. One can then imagine a network to begin with no performance 
degradation (equally useful links forming an "ordered" network), but progressively develop 
some form of disorder in which each link becomes costly to use (with cost or weight w). 
In other words, the network acquires a distribution of link weights P a (w) where a is the 
disorder parameter [9] . Once disorder has set in, new "optimal" paths (paths of least total 
weight) between each pair of nodes must be found, and a new disorder induced Limited 
Path Percolation can be defined on the basis of the path length change between the original 
unweighted network and the subsequent weighted network, reflecting a scenario more akin 
to many real networks. 

In this article, we study disorder induced LPP, and find that the model leads to a first 
order phase transition for both lattices and Erdos-Renyi networks, indicating a catastrophic 
failure of the system when the disorder exceeds the tolerance limit of the network. For Erdos- 
Renyi networks, the failure appears to be catastrophic, as our numerical results indicate 
that the phase transition occurs between an almost fully connected network to a network 
of fractional size zero in the thermodynamic limit. Through the use of theoretical results 
that explain the behavior of optimal paths, we are able to predict the location of the LPP 
transition in both network structures as a function of the disorder, and structural properties 
of the system. Our model also displays universal features with respect to the disorder, as 
demostrated by our numerical results. 

Model and Methods. 

Formulation of Limited Path Percolation induced by disorder. To be concrete, LPP in 
Ref. [5] , where link failure corresponds to link removal, is formulated in the following way: if 
a pair of nodes i and j is connected through a shortest path of length tij (number of links) 
at p = 1 (no links removed) , and of length 1^ (p) at p < 1 (fraction 1 — p of links removed) 
, i and j are considered reachable if t i j (p) < riij , where r is the tolerance factor which lies 
on the range between 1 and oo [8]. The number of reachable pairs, 5, is measured through 

S = J2^£ m -£ m ) (1) 

where o is an "origin" node chosen to minimize its impact on S (see below), and 9 is the 
Heaveside step function 9{x) = 1 is x > 0, and otherwise. S depends on both p and r, and 
the LPP phase transition occurs for the combination of these parameters at the threshold of 
the relation S ~ N. For fixed r, there is a threshold p = p c (r) at which S ~ N is achieved. 
Alternatively, given p > p c [10] . there is a r = T c (p) such that S ~ N. When r — > oo, 
p c (r -> oo) =p c . 

For disorder induced LPP, the set up is similar, but instead of considering each link to 
be kept with probability p, we change each link weight from 1 to to drawn from a random 
uncorrelated distribution P a (w) where a is the disorder parameter (defined below) [11] . The 
path lengths change from £y to lij(a) where now the later corresponds to the length of 
the optimal path between i and j. Reachability is then defined as in Ref. |S]: i and j are 
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considered reachable if (a) < riij . The number of node pairs S that remain reachable 
is a function of a, t, and N. We search for the LPP threshold by imposing S ~ N. If a 
threshold exists for a fixed tolerance t, there is a critical disorder a = a c (r) or, vice versa, 
a critical tolerance r = r c (a) that depends on the disorder parameter a. Another way to 
phrase this is to say that r and a are control parameters in disorder induced LPP, as r and 
p are in regular LPP. 

To consider a tolerance to the path length increase without considering an associated 
tolerance to the path weight increase may at first seem arbitrary, but in fact is well justified 
in that the overall weight of a path is asymptotically proportional to its length under the 
conditions of disorder we study here. Hence, choosing path length tolerance does not affect 
the qualitative nature of our results, and one choice of tolerance can be mapped onto the 
other. In concrete terms, if we imagined disorder corresponding to something like time of 
travel through a link, then total average travel time scales linearly with the trip distance. 

The general algorithm used to measure disorder induced LPP is the following. First, 
we select the ensemble of networks Q of interest. In this article we focus on 2-dimensional 
square lattices and Erdos-Renyi (ER) graphs [12] . Square lattices have no randomness, of 
course, but can be formally viewed as an ensemble with one single realization. For each 
network realization G G Q, in which all links have weight 1, we determine the path length 
Id between nodes o and i for all i ^ o. On the same network realization G (the same 
nodes and links), disorder is introduced by changing the weight of each link from I to Wij. 
Subsequently, the optimal paths between o and all other nodes i of G are determined, and 
their lengths l oi recorded. S is calculated by using Eq. (p}. To determine path lengths, we 
use the Dijkstra algorithm [13] , 

Disorder. We consider disorder distributions P a {w) characterized by a single disorder 
parameter which, for convenience, we label as a. Generally, we would expect that disor- 
der induced LPP changes as a function of the form of P a (w). However, recent work [14] 
shows that large classes of disorder distributions are essentially equivalent in the optimal 
path problem, provided a certain characteristic length scale £(a) associated with P a (w) is 
conserved (see below). In practical terms, this means that disorder distributions of different 
functional forms but with equal £(a) lead to optimal path distributions that scale in the 
same way [14]. 

This result allows us to choose a distribution that is convenient and well understood. 
Thus, we use 

P a (w) = (aw)-\ [we[l,e a ]} (2) 

for which a large amount of research has been conducted in the context of the optimal path 
problem [15H19) . Determining £ for a given distribution is addressed in |14j . and we return 
to this in the discussion of results. 

Results and Discussion. — To characterize LPP effectively, we first measure 0(S|t, a, iV),| 
the distribution of sizes of the cluster containing o, for networks of size N with tolerance 
t over disorder and network realizations. Previously [8], the LPP transition was found 
by calculating (S) = ^2 S SQ(S) and determining the parameter values at the threshold of 
(S) ~ N. Here we develop a more systematic approach, in which we look at the entire 
distribution Q(S\t, a, N) in order to explore whether disorder induced LPP exhibits a phase 
transition, and if it does, what is the order. To study the thermodynamic limit, we find it 
useful to analyze the fractional mass a = S/N, and hence ®(cr\T,a,N) (or Q(a\r,a,L) for 
lattices) . 

On a square lattice of equal sides L and N — (L + l) 2 nodes, we measure 0(<t|t, a, L) 
from node o located at the center of the lattice (x = [[(i+l)/2]], y a = [[(L + l)/2]]) with [[.]] 
indicating the next lowest integer of the argument. In Fig. [U(a), we show 0(<t|t, a — 10, L) 
with several L and t. The first interesting feature is the narrow shape of the distribution, 
indicating the presence of a characteristic mass S for given r and L. This suggests focusing 
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on the most probable value of a, labeled a* (r, a, N), i.e. 9(er* |r, a, AT) > 6(ct|t, a, A 7 ") for all 
a 7^ (7*. Also, we observe that for small r, as the system size increases, a* systematically 
decreases, but in contrast, for large r, a* increases. Between these two cases, we find a 
t for which cr*(r) remains virtually constant. For such r, labeled Tc ( Latt ' 5 S ^ N since 
cr*(r c ) = const., signaling the appearance of a phase transition. From Fig. HJa) we observe 
that cr*(ri Latt ' ) ) is close to 0.20, suggesting a first order transition. 

A systematic study of ct*(t, a, N) can be carried out with the purpose of understanding 
in more detail the phase diagram of the model. In Fig.QJb) we present a* for a = 10 and a 
range of L and r. The value of <7* was estimated from 0(ct|t, a, L) by finding an appropriate 
cubic fit for the peak of the distribution and calculating the location of its maximum. Three 
main regimes can be observed for r above, close to, and below rJ L . The determination 
of 7-c a as a function of a from simulations is done by inspection and requires exploring 
a range of r with considerable precision (St ~ 10~ 2 or even 10 -3 , becoming more sensitive 
for large a) around a certain region, in plots such as Fig. QJb) . For small r (say, close to 1 
and well below ri Latt ' 1 ), it is interesting to see that as the system size increases, a* ~ L~ 2 
indicating that the fraction of the system that is reachable is smaller than any power of N. 
The situation for r < T(? ( Latt ' j s n ot as clear: in LPP due to link removal [5] there is a regime 
of power law sizes of S, whereas here such regime is not evident for lattices but seems to be 
present for ER networks (see below); our current theory does not shed light on the matter. 
For r > Tc ( Latt ) 5 we fi nc [ a progressive increase of a* with respect to L, with a saturation 
value that depends on t; a* gradually approaches 1 as r —> oo. 

To measure 0(<t|t, a, N) in ER networks, we sample over network realizations of G £ Q, 
and for each G choose an origin o at random. In Fig. [T](c) we present the relevant simulation 
results. The qualitative features of 0(er|r, a, N) for lattices are also present in ER networks, 
including the existence of a critical t, labeled ri ER '. In contrast to lattices, a*(Tc E ^) is 
close to 1, indicating a very dramatic LPP transition, in which a slight change of t around 
ri ER ' 1 leads to a transition between almost entirely reachable to entirely unreachable global 
network states. We also observe that for r < ri ER - ) , there seems to be a power-law decaying 
relation between a* and N, consistent with a fractal size object below the threshold, with 
the decay exponent r-dependent. 

To analyze the problem further, we define the quantity ay = ^j/^ij, called the length 
factor, for each node pair ij [20], which measures the fractional increase of the path between 
i and j, and explore the distribution $(a) and its cumulative F(a) — J $>(a )da over 
realizations of P a (w) (and Q for ER networks). Note that the LPP reachability condition is 
ay < t. 

Figure HJa) shows the measurement of F(a\a,L) in lattices, and in the inset the distri- 
bution <$>(a\a,L), which is a well concentrated function around its maximum a — a c . The 
cumulative F(a) increases sharply around a = a c rapidly approaching 1, which indicates 
that many node pairs satisfy a < a c . The sharpness of $(a) increases for larger L, and F(a) 
becomes more step like, while a c remains in the same location. Note that a c can also be de- 
termined at the location where F(a\a,L) for increasing L cross over each other. Increasing 
a, on the other hand, leads to increasing a c , consistent with path lengths becoming longer 
under more disordered conditions. In Fig. [He), we focus on F(a\a, N) for ER networks, 
and observe similar features to those in the case of lattices (we display the cumulative only 
as the small values of path lengths in ER graphs produce large discretization fluctuations), 
apart from the asymptotics which appear to be slower. 

Scaling of path length as a function of disorder. In order to understand the previous 
results, we consider the current knowledge on the problem of optimal paths, which has 
received considerable attention in the context of surface growth and domain walls [T4 [IT5f2T] 
in the physics literature. For the purpose of clarity we briefly review these results here, 
starting with lattices and extending the discussion to networks. 
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Disordered lattices often exhibit optimal paths which are self-affine, characterized by 
lengths which scale linearly with , with a constant prefactor dependent on the roughness 
exponent related to the disorder [3T]. This limit of self-affine paths has become known as 
the weak disorder limit. 

Another scaling regime has been recognized |15j when the disorder approaches the so- 
called strong disorder limit. In this limit, each link weight in the network is very different 
to any other link weight, progressively forcing the optimal paths to lie inside the minimum 
spanning tree where their lengths scale as l£ px , d op t being the scaling exponent of the 
shortest path in the minimum spanning tree |14[|16l[2"2"] . 

A general theory explaining these optimal path limits points out that weak and strong 
disorder are separated by a disorder length-scale £ which depends on the disorder distribution 
P a (w) and some lattice dependent features [UJ[T7]- Optimal paths covering a distance 
smaller than £ are in strong disorder, and those covering a larger distance are in weak 
disorder (provided the system is large enough so that ( « L). Thus, £ is the weak-strong 
disorder crossover length. The weak and strong disorder scaling regimes for I- can be 
expressed by the scaling relation 

/ ^^ff^A f( T \-f x > £>l(weak) 
« ? Vw' \ const., x <C I (strong). {) 

The length-scale £ can be determined by the relation £ = [p c /(w c P a (w c ))] L ' [T3], where p c 
is the percolation threshold of the lattice, v the correlation length exponent of percolation, 
and w c is the solution to the equation p c = J P a (w)dw, i.e., the weight for which the 
cumulative distribution of P a (w) is equal to p c - 

Based on the previous arguments, we can now postulate the properties of Q(a). We 
concentrate on the weak disorder limit because it is the only possible regime in which 
an LPP transition could take place 23 . In this regime, based on Eq. ©, we find that 

~ £ dopt (iij l '£) = £ dopt_1 %, where the parenthesis corresponds to £jj in the scale of the 
crossover length £, and £ dopt to the length of the path within this crossover length (see 
Fig. 3 in Ref. [17]). This relation indicates that the typical a is given by a c ~ ^ij/^ij = 
£ > d ° Ft ~ 1 lij / lij — ^ opt_1 . For P a (w) of Eq. (jSJ, £ = (ap c ) v calculated according to the above 
formalism, producing 

«i Latt) ~ (ape)"****- 13 . (4) 

In order to test this, we present in Fig.[3ja) (main) the scaled curves &(a/a c ) for various L 
and a, where a c is taken from Eq. ((4]). The collapse is consistent with our scaling picture, 
and indicates that indeed there is a clear path length increase (ap c ) l/ ' d ° pt-1 - ) that explains 
the empirical results. 

Given the large fraction of node pairs for which a is close to a c , we postulate that 

Tc (Latt) = a (Latt) (fl) _ (op,.)^***- 1 ), (5) 

i.e., the tolerance necessary to obtain the LPP transition is equivalent to the most probable 
length factor. To test this relation, we find by inspection the values of r ^ Latt ) as a function of 
a, and plot them in Fig.[3ja) (inset). The relation between Tc ( Latt ) anc [ a can b e n t to a power 
law with exponent 0.297±8, which is very close to the predicted v(d opt — 1) = 0.293, strongly 
supporting Eq. ([5]). In addition, we find that a* (r c (Latt) ) - 0.2 ±0.04, independent of a. The 
previous results indicate that a* is universal. We have tested this by comparing the shape of 
ct*(t, a, L) for different a, r, and L, and have found that adjusting for a given combination 
of these values, the curves for a* can be made to overlap, supporting universality. 

For ER networks, it is known that their path length distributions are concentrated due 
to their random structure, leading to a large number of lengths being similar to an overall 
typical length [IB] . Thus, we simplify our analysis by focusing on the typical lengths before 
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and after the introduction of disorder. Before disorder sets in, the typical path length is 

Per ~ ; — 7TT- (6 
log(fc) 

However, after weak disorder sets in, it becomes 

?-ER ~ V a Vc (7) 

V a Pc J 

which emerges from the relation that is the equivalent to Eq. ([3]) now applied to networks. 
fi is a quantity not yet characterized in the literature, and depends on (k) and also weakly 
on ./V and a; for our simulations, fi « 4. The ./V-dependent ratio between £er and t ER , 
which we label /3(N) (slightly different to a because the later applies to each pair of nodes, 
but the former to the overall typical distances) , is given by 



1 log(ap c ) 



3 log TV 



a^-e^(N) (8) 



where ai™^ = p,ap c log(fc) /3, the asymptotic value of /?, and e( ER ) (N) = \iap c log(fc) log(ap c )/ log 
a finite size correction that vanishes logarithmically with N as N — » oo. In analogy to lat- 
tices, we hypothesize that Tc ER \n) scales as /3(N), which includes the size corrections with 
N. Also, in the limit N —> oo, we define ri™" 1 = Oc^» • 

To test Eq. ([8j), we introduce the rescaled variable a'^ — (aij+e( ER \N))/oi^\ and plot 
F(a \a, N) in Fig. |3jb). The collapse of the curves is excellent, and supports the validity of 
our assumptions. In the thermodynamic limit we expect 

(ER) = (ER) „ MQPclQg(fc) 

c,oo W C,00 5 W ) 

but it is important to keep in mind the large finite size corrections to ri ER ' ) , which make it 
more similar in value to Eq. (|8|) for finite N. 

Conclusions. — We study disorder induced Limited Path Percolation and determine 
that a percolation transition occurs for a critical tolerance r c . The critical tolerance increases 
together with the heterogeneity of the disorder. Numerical results indicate that LPP displays 
universality. Also, the phase transition is first order, with the discontinuity in the order 
parameter cr*(r c ) independent of the disorder. For lattices cr*(r c ) is of the order of 0.2; 
for Erdos-Renyi networks, our numerical results suggest it may approach 1 signaling a 
catastrophic transition. The concept of length factor applied to the theory of optimal 
paths predicts a typical factor a c which, in turn, predicts r c . We believe the tolerance 
thresholds predicted here reflect reachability conditions of some real networks under real 
failure scenarios (congestion, maintenance, etc.) better than regular percolation. 
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Fig. 1: (a)0(a|r, a = 10, L) vs. a for the square lattice for various r (in legend) and L = 160, 320, 640 
(increasing symbol size corresponds to increasing system size). Note that for r = 1.4 the peaks of 
the distributions a* remain almost in the same location with increasing TV indicating that r j LaU ' 
is near this value. For r < ri Latt ', a* decreases with TV, and for r > r c ' Latt ', a* increases with TV. 
Further detailed simulations reveal a better estimate of r ^ Latt ' [ s i.3g for a = 10. From Q(a\r, a, L), 
we determine a* by cubic regression of the top points of the peak of O. (b) a* (r, a = 10, L) vs. L for 
square lattice, for several r and L. With r = T( ( Latt ) _ i 3g ^ e curv e stays constant indicating the 
phase transition; for r approaching 1 (path lengths cannot increase much), a* scales as L~ 2 (solid 
line), indicating that the largest reachable component is logarithmic in size, (c) Q(a\r, a = 8, TV) 
vs. a for Erdos-Renyi networks of (k) = 3, TV = 11000,59000,307000 (increasing symbol size 
corresponds to increasing system size), and various r. In this case, for r ~ 3.7, a* stays roughly 
constant with increasing TV, indicating rj ER ' is close to this value. A more detailed analysis reveals 
that ri ER ' = 3.74 for a — 8 and (k) — 3. (d) <t*(t, a = 8, TV) for ER networks, with r values 
indicated in the legend, and TV ranging from 11000 to 307000. It is clear that a* increases with TV 
for t > ri ER ' and decreases for r < i~i ER ' ) . 
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Fig. 2: Cumulative F(a\a, L) vs. a for 2-dimensional lattices, and F(a\a, N) for ER networks. Plot 
(a) for lattices corresponds to different L and a. The inset shows $(a|a, L) with logarithmic scale in 
the vertical axis, emphasizing the fast decay of the distribution which is even more pronounced for 
larger system sizes. This supports the idea that, at the asymptotic limit L — > oo, the distribution 
is highly concentrated around a c (highlighted by the arrow). Note that the location of a c also 
corresponds to the location where F(a\a, L) cross over each other as L increases. For increasing a, 
the location of a c shifts to the right. Panel (b) for ER networks corresponds to F(a\a, N) for several 
N and a (distribution $ is omitted due to discretization fluctuations). a c is such that F(a\a,N) is 
approaching 1. Increasing a leads to distributions shifted to the right. The inset shows the crossing 
over between distributions with fixed a and increasing N. 
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Fig. 3: Testing Eqs. ([5| and ([8|. (a) <E>(a/a c ) and F(a/a c ) vs a/a c for lattices with a c ~ 
(opc)"^^^ 1 ', for various a and L specified in the legend. Both $ and F collapse supporting 
the hypothesis that a c responds to Eq. Q for lattices. Inset: r j Latt ) VSp a from data and least 
squares fit. The relation predicted by our theory is r c ~ (ap c )"' dopt_1 ' which is equal to 0.293 from 
independent work, within error of 0.297 ± 0.008 obtained from the fit from our simulation results, 
(b) F(a \a,N) vs. a for ER networks for a combination of and a. The collapse of the curves 
strongly supports Eq. ((HJ> . 
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